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Abstract 

The  heat  and  species  transport  processes  in  a  tubular  type  solid  oxide  fuel  cell  (SOFC)  that  works  in  a  cell  stack  were  analyzed  and 
modeled.  Since  most  of  the  single  tubular  SOFCs  working  in  a  cell  stack  share  the  same/similar  chemical/electrochemical  and  heat/mass 
transfer  conditions,  it  is  plausible  to  assume  that  heat  and  species  are  not  exchanged  between  one  cell  and  its  neighboring  cells.  Therefore, 
a  surrounding  fuel  flow  space  was  outlined  controllable  by  a  specific  single  cell,  for  which  zero  flux  was  assumed  at  its  boundary  in 
neighborhood  with  other  cells.  The  numerical  model  subjects  such  a  cell  and  its  controllable  fuel  flow  space  to  a  two-dimensional  analysis 
for  the  flow,  heat/mass  transfer  and  chemical/electrochemical  performance.  Computations  were  performed  for  three  different  tubular  SOFCs 
having  practical  operating  results  available  from  publications  by  different  researchers.  The  numerical  results  of  the  terminal  voltages  for 
those  different  SOFCs  showed  very  good  agreement  with  the  published  experimental  data.  It  is  expectable  that  the  proposed  numerical 
model  be  used  to  significantly  help  the  design  and  operation  of  a  SOFC  stack  in  practical  applications. 
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1.  Introduction 

A  solid  oxide  fuel  cell  (SOFC)  is  one  of  the  most  promis¬ 
ing  type  of  fuel  cells  currently  being  considered  to  be  used  as 
a  power  source  for  automobiles  and  stationary  power  plants 
[1,2].  A  SOFC  has  a  oxide-ion-conductive  electrolyte  and 
works  by  the  oxidation  of  a  fuel  at  the  anode  side.  Both  car¬ 
bon  monoxide  and  hydrogen  can  serve  as  fuels,  which  pro¬ 
vides  the  SOFC  with  an  advantage  in  fuel  adaptability.  For 
maintaining  the  high  oxide-ion  conductivity  of  the  solid  ox¬ 
ide  electrolyte,  SOFCs  must  operate  at  a  high  temperature 
of  800-1000  °C.  This  also  provides  SOFCs  with  another  ad¬ 
vantage  in  that  natural  gas  can  be  reformed  in  the  cell  stack 
[3,4]. 

A  SOFC  essentially  consists  of  two  porous  electrodes 
separated  by  a  dense,  gas-tight,  oxide-ion-conducting  elec¬ 
trolyte.  Planar  and  tubular  geometries  are  the  two  popular 
designs  for  solid  oxide  fuel  cells.  The  technique  of  designing 
and  manufacturing  a  tubular  type  SOFC  is  relatively  mature, 
and  it  is  possible  to  construct  a  medium  and  large  scale  cell 
stack  for  automobile  and  stationary  power  generation.  There 
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have  been  growing  experimental  test  results  on  the  overall 
electrochemical  performance  for  tubular  type  SOFCs.  How¬ 
ever,  it  is  still  far  from  enough  to  guide  the  design  and  oper¬ 
ation  of  a  solid  oxide  fuel  cell,  especially  the  cell  stack.  To 
meet  the  need  of  understanding  the  operational  details  and 
the  internal  temperature  and  species’  mass  fraction  parame¬ 
ters  when  a  fuel  cell  works  in  a  cell  stack,  the  present  study 
develops  a  reliable  numerical  model  validated  by  available 
experimental  data.  Localized  experimental  measurements  in 
a  SOFC  working  at  high  temperature  is  difficult,  so  numer¬ 
ical  analysis  is  the  most  viable  way,  at  the  present  stage,  for 
studying  the  inside  details  of  heat  sources,  mass  transport 
and  temperature  distributions. 

There  are  some  valuable  studies  which  model  the  tubular 
type  solid  oxide  fuel  cells  [5-9].  One  obvious  weakness 
of  these  studies  is  that  they  did  not  make  any  analysis  for 
the  SOFC  performance  when  it  is  working  in  a  cell  stack. 
Although  it  is  likely  that  a  SOFC  might  demonstrate  similar 
performance  whether  it  is  operated  as  a  single  SOFC  unit  or 
in  a  cell  stack,  a  pertinent  modeling  or  reliable  analysis  for 
this  is  very  important  for  the  practical  operation  of  a  SOFC 
stack.  The  primary  effort  in  this  study  will  be  focused  on 
this  issue. 

Another  weakness  of  the  previous  modeling  work  by 
other  researchers  is  that,  instead  of  employing  a  field  so- 
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Nomenclature 

^  cell 

area  of  electrolyte  layer  (m2) 

Cp 

heat  capacity  (J/(kg°C)) 

D\,m 

mass  diffusivity  of  gas  I  in  gas  mixture 
(m2/s) 

E 

electromotive  force  (V) 

F 

Faraday’s  constant  96486.7  (C/mol) 

AG 

Gibbs  free  energy  change  (J/mol) 

A  G° 

standard  Gibbs  free  energy  change  (J/mol) 

H 

the  interval  from  the  exit  of  air-inducing 
tube  to  the  closed  end  of  cell  (m) 

AH 

enthalpy  change  of  the  reaction  (J/mol) 

i 

current  (A) 

^cell 

average  current  density  based  on  area  of 
electrolyte  layer  (mA/m2) 

L 

length  of  fuel  cell  (m) 

m 

mass  flux  (g/(m2s)) 

M 

formula  weight  (g/mol) 

Mf 

mole  flow  rate  of  fuel  (mol/s) 

P 

pressure  (Pa) 

P 

ratio  of  a  pressure  over  1.013  x  105  Pa; 
center  position  of  a  control  volume  or  area 

k 

volumetric  heat  source  (W/m3) 

Q 

heat  generation  (W) 

r 

radial  coordinate  (m) 

ra,  rc,  re 

average  radius  of  the  anode,  cathode 
and  electrolyte,  respectively  (m) 

R 

universal  gas  constant  8.31434  (J/(molK)) 

Ra,  Rc 

discretized  electrical  resistance  in  the 
anode  and  cathode,  respectively  (£2) 

Re 

discretized  ionic  resistance  in 
electrolyte  layer  (Q) 

AS 

entropy  change  of  the  reaction 
(J/(mol  K)) 

T 

temperature  (K) 

u ,  V 

velocities  in  axial  and  radial  directions, 
respectively  (m/s) 

U 

utilization  factor  for  hydrogen  and 
oxygen  (%) 

V1, 

electrical  potential  at  anode  and 
cathode,  respectively  (V) 

Tcell 

cell  terminal  voltage  (V) 

coordinate  in  axial  direction  (m) 

Ax 

mesh  size  or  a  section  in  axial 
direction  (m) 

X 

mole  fraction  (%) 

y,  z 

mole  numbers  of  CH4,  CO  and  H2  that 
involve  in  the  reforming,  shift  and  electro¬ 
chemical  reactions,  respectively  (mol/s) 

Y 

mass  fraction  (%) 

Greek  letters 

<5a,  <5C,  <5e 

thickness  of  anode,  cathode  and 
electrolyte  layers,  respectively  (m) 

^act 

activation  polarization  (V) 

AO 

mesh  size  in  peripheral  direction  (°) 

A 

thermal  conductivity  (W/(m°C)) 

/X 

dynamic  viscosity  (kg/(m  s)) 

P 

density  (kg/m3) 

PQ 

resistivity  (£2  m) 

Subscripts 

air 

air  flow 

ch4 

methane 

CO 

carbon  monoxide 

C02 

carbon  dioxide 

e,  w,  n,  s 

interfaces  between  the  P-centered 
control- volume  and  its  neighboring 
control- volumes  of  east,  west,  north 
and  south 

E,  W,  N,  S 

control- volumes  that  locate  east,  west, 
north  and  south  to  the  P-centered 

control-volume 

f,  fuel 

fuel 

h2 

hydrogen 

h2o 

water  vapor 

I 

gas  species 

02 

oxygen 

P 

center  of  a  control-volume  or  -area 

PR 

equilibrium  of  reforming  reaction 

PS 

equilibrium  of  shift  reaction 

V 

axial  position 

Ax 

axial  section 

Superscripts 

a 

anode 

air 

air  flow 

c 

cathode 

ch4 

methane 

CO 

carbon  monoxide 

co2 

carbon  dioxide 

e 

electrolyte 

f,  fuel 

fuel 

h2 

hydrogen 

in 

inlet  of  fuel  or  air  in  a  concerned  section 

I 

gas  species 

02 

oxygen 

out 

outlet  of  fuel  or  air  in  a  concerned  section 

R 

reaction 

lution  for  the  flow,  heat  and  mass  transfer  for  obtaining 
the  local  species’  mass  fractions,  temperatures  and  the 
electromotive  forces  in  the  fuel  cell,  most  authors  used  con¬ 
stant  heat  and  mass  transfer  coefficients  based  on  a  fully 
developed  laminar  flow  approximation  at  constant  wall 
temperature  and  mass  flux.  In  fact,  the  tubular  SOFC  is  a 
heat- generating  tube  with  different  flow  streams  on  both  the 
inner  and  outer  side.  The  coupled  heat  and  mass  transfer  in 
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the  two  streams  are  not  applicable  to  the  case  of  constant 
wall  temperature/concentration  or  the  case  of  constant  wall 
flux  of  heat/mass.  Therefore,  the  fully  developed  laminar 
flow  approximation  for  Nusselt  and  Sherwood  numbers  at 
constant  wall  temperature/concentration  can  no  longer  be 
applied  to  the  SOFC.  There  is  not  doubt  that  significant  de¬ 
viation  might  occur  as  has  been  well  discussed  in  [10].  As 
has  been  discussed  in  the  author’s  previous  work  [11],  we 
will  avoid  the  fully  developed  laminar  flow  approximation 
of  Nusselt  and  Sherwood  numbers  at  constant  wall  temper¬ 
ature/concentration  in  this  study,  instead,  a  complete  field 
solution  of  the  governing  equations  for  the  heat  and  mass 
transfer  in  the  entire  domain  of  a  SOFC  working  in  a  cell 
stack  will  be  employed. 


Simplified 
cross-sectional 
domain  of  one  cell 


(0)  Cathode 


Air  in  ^ 

Interconnection 


2.  Model  development 


Fig.  2.  Orientation  of  a  tubular  SOFC  in  a  cell  stack. 


2.1.  Delineation  of  the  domain  pertain  to  a  single  SOFC 
in  a  cell  stack 

Fig.  1  shows  the  schematic  view  of  a  typical  tubular 
SOFC.  The  layers  of  cathode,  electrolyte  and  anode  are  lam¬ 
inated  onto  the  gas  diffusible  porous  support  tube.  The  cell 
tube  is  fabricated  in  the  structure  like  a  test  tube  so  that 
air  is  supplied  through  a  concentric  air-inducing  tube  in¬ 
serted  inside  the  tubular  cell  from  its  open-end.  Oxygen  in 
air  is  ionized  at  the  cathode  (air  electrode),  and  the  oxide 
ions  produced  are  conducted  through  the  electrolyte  to  the 
anode  (fuel  electrode).  Oxidization  of  the  fuel  (hydrogen 
or  carbon  monoxide)  that  diffuses  from  the  core  region  of 
the  fuel  stream  occurs  at  the  anode.  In  the  meantime,  the 
products  from  the  chemical/electrochemical  reactions  dif¬ 
fuse  back  into  the  fuel  stream  and  are  removed  together  with 
the  remaining  unreacted  fuel. 

In  a  practical  tubular  SOFC  stack,  many  tubular  cells  are 
mounted  in  a  container  in  a  cell  bundle,  as  shown  in  Fig.  2. 
For  most  of  the  tubular  cells  mounted  in  the  bundle,  each  of 
them  is  surrounded  by  four  others.  Therefore,  it  is  possible, 
and  very  likely,  that  most  of  the  single  tubular  SOFCs  work 
under  the  same  environment  of  temperature  and  concentra¬ 
tions  of  gas  species.  This  allows  us  to  define  a  controllable 


domain  in  the  cross-section,  which  pertains  to  one  particular 
single  cell,  as  outlined  by  the  dashed-line  square  in  Fig.  2. 
There  must  be  no  flow  velocity  and  fluxes  of  heat  and  mass 
across  the  outline.  This  will  significantly  simplify  the  anal¬ 
ysis  of  a  cell  stack.  In  other  words,  through  analysis  of  the 
heat/mass  transfer  and  the  chemical/electrochemical  perfor¬ 
mance  of  one  single  cell  and  its  controllable  area,  one  can 
obtain  results  that  are  very  useful  for  evaluating  the  perfor¬ 
mance  of  a  cell  stack. 

With  the  longitudinal  direction  also  being  considered,  the 
heat  and  mass  transfer  in  the  above  outlined  square  area  en¬ 
closing  the  tubular  SOFC  are  in  three  dimensions.  For  a  so¬ 
lution  of  the  three-dimensional  governing  equations  of  mo¬ 
mentum,  energy  and  species  conservation,  a  large  number  of 
discretized  mesh  is  necessary  and  that  will  make  the  compu¬ 
tation  too  difficult  to  be  conducted  with  a  personal  computer. 
In  order  to  reduce  computation  cost,  the  square  area  enclos¬ 
ing  the  tubular  SOFC  is  approximated  to  be  an  equivalent 
circular  area  and  therefore,  the  domain  enclosing  the  inter¬ 
ested  single  tubular  SOFC  is  viewed  as  a  two-dimensional 
axial- symmetric  one,  as  seen  in  Fig.  2.  However,  the  zero 
flux  or  insulation  of  heat  and  mass  transfer  at  the  boundary 
remains  unchanged  even  though  this  geometric  approxima¬ 
tion  is  applied. 


Fig.  1.  Schematic  of  a  tubular  solid  oxide  fuel  cell  (SOFC). 


2.2.  Electro  chemical/ chemical  reactions  and 
electromotive  force 

If  the  SOFC  operates  based  on  pre-reformed  natural  gas 
[12],  the  fuel  stream  might  have  components  of  CH4,  CO, 
CO2,  H2,  and  H2O.  Then,  the  reforming  and  shift  reactions 
will  occur  along  with  the  electrochemical  reaction,  i.e.  the 
oxidization  of  hydrogen 

CH4  +  H2O  CO  +  3H2  (reforming)  (1) 

CO  +  H20  C02  +  H2  (shift)  (2) 

h2  +  \o2  H2O  (electrochemical  oxidation)  (3) 
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In  fact,  the  oxidization  of  CO  also  occurs  at  the  anode. 
However,  the  variation  of  CH4,  CO,  CO2,  H2  and  H2O 
must  satisfy  the  chemical  reaction  equilibrium  of  Eqs.  (1) 
and  (2),  therefore,  there  is  no  necessity  to  distinguish 
whether  the  electrochemical  reaction  step  involves  oxida¬ 
tion  of  H2  or  CO  because  both  the  oxidation  of  H2  and  CO 
yield  the  same  electromotive  force  [2]. 


-AG  _  —A G°  RT  /  Ph2  Pf  \ 

2  F  =  2  F  +2Fln  \  PH2o  ) 


where  A G°  is  the  variation  of  standard  state  (1 .013  x  105  Pa) 
Gibbs’  free  energy  of  the  reaction  of  Eq.  (3),  Pu2  and  Ph2o 
are  the  ratios  of  partial  pressure  over  the  standard  state  pres¬ 
sure  of  1.013  x  105  Pa  for  hydrogen  and  water  vapor  at  the 
electrolyte/anode  interface,  respectively,  and  Po2  is  that  for 
oxygen  at  the  electrolyte/cathode  interface. 

In  the  electrochemical  reaction  proceeding  at  temperature, 
T,  the  changes  of  chemical  enthalpy,  AH ,  entropy,  AS  and 
Gibb’s  free  energy,  AG,  have  the  thermodynamic  relation¬ 
ship: 

AG  =  AH  -  TAS  (5) 


Theoretically,  —AG  is  converted  into  electrical  power  as 
was  given  in  Eq.  (4),  while  —TAS  is  released  in  the  form  of 
thermal  energy  during  the  electrochemical  reaction.  Along 
with  the  consumption  and  production  of  the  gas  species  in 
fuel  and  air  flows,  the  mass  fractions  of  the  reactants  and 
products  from  the  electrochemical/chemical  reactions  also 
vary  in  the  flow  streams  of  fuel  and  air.  This  yields  the 
localized  electromotive  force,  E ,  over  the  electrolyte  layer, 
which  yields  the  localized  ion/charge  transfer  rate  through 
the  electrolyte  layer 


E  -  rfci  ~  (Vc 

Rt 


where  Vs  and  V0  are  the  potentials  in  anode  and  cathode, 
respectively,  ?7act  is  the  over-potential  incurred  by  the  ac¬ 
tivation  polarization  [13]  which  occurs  at  both  anode  and 
cathode, 

2 RT  /  /  \ 

>7  =  — psmh  I  -=-  )  (7) 

neF  \2ioJ 


where  ne  is  electrons  transferred  per  reaction  (for  the  re¬ 
action  of  Eq.  (3),  it  is  2),  i  is  the  current  density  across 
electrolyte  and,  /o  is  the  exchange  current  density  which  is 
5300  A/m2  for  anode  and  2000  A/m2  for  cathode.  The  Re  in 
Eq.  (6)  is  the  ionic  resistance  of  the  electrolyte  layer  in  the 
thickness  of  8e  and  a  unit  area  of  A  A,  which  is  in  the  form  of 


where  the  yOe  is  the  ionic  resistivity  of  electrolyte,  which  is  a 
strong  function  of  temperature  [14,15]  as  given  in  Table  1. 

In  a  tubular  type  SOFC,  the  current  is  circumferentially 
collected  to  the  nickel  felt.  The  ohmic  loss  in  the  circumfer¬ 
ential  pathway  is  significant  and  certainly  not  negligible.  In 


Table  1 


Properties  of  fuel  cell  components 


Thermal 
conductivity 
(W/(m  K)) 

Electric  or  ionic 
resistivity  (£2  cm) 

Cathode 

6.0 

0.0081 14  exp(500/7) 

Electrolyte 

2.7 

10.0exp[10092(1.0/r  -  1.0/1273.15)] 

Anode 

11.0 

0.00298  exp(-l  392/7) 

Support  tube 

1.1 

— 

Air-inducing 

1.1 

— 

tube 

this  work,  we  apply  a  network  circuit,  as  shown  in  Fig.  3,  to 
the  three  layers  of  anode,  electrolyte  and  cathode  to  analyze 
the  current,  or  charge  transfer  rate.  Because  the  current  col¬ 
lection  is  symmetric  in  the  peripheral  direction  in  the  cell 
tube,  only  one-half  tube  is  deployed  (over  14  nodes  in  pe¬ 
ripheral  direction)  in  the  analysis.  Applying  Kirchhoff’s  law 
of  current,  we  obtain  the  equation  associating  the  potentials 
of  the  interested  grid  P  with  the  potentials  of  its  neighboring 
points,  east,  west,  north,  south  and  the  corresponding  grid  P 
in  the  cathode 


ya  _  ya  va  —  ya 
VE  v?  vw  v? 


R* 

iVe 


R* 

jfVw 


ya  _  ya  va  —  ya 
_l_  I  VN  VP  /  S  V? 


+ 


-  ypa  -  (£P  -  >7act) n 


R* 

j£Vn 


=  0 


R* 


In  the  same  way  for  a  grid  P  in  the  cathode,  there  is 


where  7?a  and  Rc  are  the  discretized  resistance  in  anode  and 
cathode,  respectively,  which  are  determined  according  to  the 
resistivity,  the  length  of  current  path  and  the  area  the  current 
acted  on.  The  material  resistivities  of  both  cathode  and  anode 
are  also  given  in  Table  1,  as  the  function  of  temperature. 

With  all  the  equations  for  the  discretized  grids  in  both 
cathode  and  anode  enclosed,  a  matrix  represented  by  the  pair 
of  Eqs.  (9)  and  (10)  is  formatted.  For  the  solution  of  such 


Fig.  3.  Allocation  of  nodes  for  computation  of  electric  potentials  in 
electrodes. 
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a  matrix  equation  for  potentials,  the  following  approxima¬ 
tions  are  useful:  (1)  at  the  two  ends  of  the  cell  tube,  there 
is  no  longitudinal  current  and,  therefore,  insulation  condi¬ 
tion  is  applicable;  (2)  at  the  symmetric  plane  A-A  as  shown 
in  Figs.  1  and  3,  there  is  no  peripheral  current  in  the  cath¬ 
ode  or  anode,  unless  the  cathode  or  anode  is  in  contact  with 
nickel  felt,  through  which  the  current  flows  in  or  out;  (3) 
the  potentials  on  the  nickel  felts  are  assumed  uniform  due 
to  their  high  electrical  conductivities;  (4)  since  the  potential 
difference  between  the  two  nickel  felts  is  the  cell  terminal 
voltage,  the  potential  at  the  nickel  felt  in  contact  to  the  an¬ 
ode  layer  can  be  assumed  to  be  zero,  thus  the  potential  at 
the  nickel  felt  in  contact  with  the  cathode  will  be  the  termi¬ 
nal  voltage  of  the  fuel  cell.  Once  we  obtained  all  the  local 
electromotive  forces  from  Eq.  (4),  the  only  unknown  condi¬ 
tion  for  the  matrix  equation  is  either  the  total  current  flowing 
out  from  the  cell  or  the  potential  at  the  nickel  felt  in  contact 
with  the  cathode.  This  provides  two  approaches  for  the  nu¬ 
merical  simulation  of  the  fuel  cell.  In  case  the  total  current 
taken  from  the  cell  is  prescribed  as  the  initial  condition,  the 
terminal  voltage  will  be  predicted  as  the  simulation  result. 
On  the  other  hand,  one  can  prescribe  the  terminal  voltage 
and  target  the  simulation  computation  in  obtaining  the  total 
current,  i.e.  the  summation  of  local  current  i  across  the  en¬ 
tire  electrolyte  layer.  In  the  present  work,  the  first  approach 
was  adopted. 

With  the  potentials  being  obtained  in  an  electrode  layer, 
the  Joule  heating  in  that  electrode  in  an  interested  area,  con¬ 
trolled  by  P,  will  be 


The  thermodynamic  heat  generation  occurring  at  the  an¬ 
ode/electrolyte  interface  in  the  interested  area,  controlled  by 
P,  is: 


(AH  —  AG)/P 
IF 


(16) 


The  current  ip  and  heating  rates  distribute  non-uniformly 
along  the  axial  and  peripheral  directions.  Because  two- 
dimensional  flow  and  heat/mass  transfer  are  assumed,  the 
peripheral  average  of  the  local  heating  rates  and  ionic  con¬ 
ducting  rate  or  current  across  the  electrolyte  is  introduced 
into  the  axial- symmetric  equations  governing  the  conserva¬ 
tion  of  momentum,  heat  and  mass. 

If  the  SOFC  operates  based  on  pre-reformed  natural  gas, 
the  variation  of  CH4,  CO,  CO2,  H2,  and  H2O  associate  each 
other  and  need  to  be  determined  based  on  the  condition  of 
chemical  reaction  equilibrium.  For  convenience  of  analysis, 
the  mole  numbers  of  CH4,  CO  and  H2  reacted  in  reactions 
by  Eqs.  (1)— (3)  in  a  concerned  section  in  fuel  channel  are 
represented  by  x ,  y  and  z,  respectively.  Then,  the  variation 
of  mole  numbers  of  all  the  species  in  an  interested  section 
of  fuel  channel  from  its  inlet  to  outlet  is: 


CHf1  =  CH41  -  3c 

(17) 

C0out  =  COin  +  3c  -  y 

(18) 

COT  =  CO"’  +  v 

(19) 

H°ut  =  nf  +  35c  +  y-z 

(20) 

+ 


H20out  =  H2Oin  -  x  -  y  +  z 


(21) 


The  variation  of  the  total  mole  numbers  of  the  fuel  flow  from 


+ 


(AxPra  A#p<5a)-1 


(11) 


inlet  to  outlet  of  the  interested  section  becomes: 

Mfout  =  Mf  +  23c  (22) 


+ 


(AxprcA#pi5c)  * 


(12) 


(E P  -  ijf  -  V£  +  V*)2 


(AxPre  A#p<5e)-1 


(13) 


where  the  r  and  8  with  corresponding  superscripts  a,  c  and 
e  are  the  average  radius  and  thickness,  respectively,  for  the 
anode,  cathode  and  electrolyte,  and  Axp  and  AOp  are  the 
P-controlled  mesh  size  in  axial  and  peripheral  directions  as 
shown  in  Fig.  3.  The  heating  rate  induced  from  activation 
polarization  in  anode  and  cathode  is  in  the  form  of: 


•  act,  a 
9[p 


act,  a 

Wp 


AxprR  A6p8R 


(14) 


act,c 

iprjp 


According  to  the  chemical  equilibrium,  the  equilibrium  con¬ 
stants  of  the  reaction  by  Eqs.  (1)  and  (2)  are  in  the  forms  of 


Pl  Pco 

Kpp  =  - =  exp 

Pch4  Ph2o 


^  ^reforming 

RT 


) 


(23) 


A'ps 


PCQ2  /J"2  =  l  AGs°hift 

PcoPr2o  \  PP 


(24) 


where  Ph2o>  Fh2,  Pch4>  Pco  and  Pco2  are  the  ratios  of 
the  partial  pressure  over  the  standard  state  pressure  of 
1.013  x  105  Pa  for  the  corresponding  species  in  fuel  stream 
of  the  interested  section.  When  expressing  the  partial  pres¬ 
sures  in  terms  of  the  species’  corresponding  mole  numbers 
and  the  overall  pressure,  the  equilibrium  constants  is  in  the 
form  of 


[(COin  -\- x  —  y)/ (Mjn  +  2x)] 

X  |  (H"’  +  31-  +  V  -  z)/(M'P  +  2x)]3  P2 

[(CHjj1  -  5c) /(Mf  +  2x)\ 

x [(H2Oin  -x-y  +  z)/(Mf  +  25c)] 


(25) 


Axprc  A0p8c 


(15) 
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[(Hi,n  +  3v  +  y  -  z)/(Mf  +  23c)] 
x  [(CO1”  +  y)/{Mf  +  2x)\ 

KPS  =  - — ^ ? -  i - - -  (26) 

[(COin  +  X  -  y)/(Min  +  23c)] 

x  [(H2Oin  -x-y  +  z)/(Min  +  23c)] 


where  P  is  the  ratio  of  overall  pressure  over  the  standard 
state  pressure  of  1.013  x  105  Pa  for  a  particular  section  in 
the  fuel  stream. 

In  a  particular  section,  the  reacted  H2  and  02  are  related 
to  the  current  transfer  rate  across  the  electrolyte  layer  of 
the  section,  which  is  obtainable  through  the  solution  of  the 
network  circuit  as  discussed  in  the  preceding  paragraphs. 


z  = 


i 

IF 


(27) 


i 

4 F 


(28) 


With  the  2  and  the  temperature  of  gas  species  available, 
the  x  and  y  are  obtainable  from  Eqs.  (25)  and  (26).  The 
reforming  and  shift  reactions  take  place  at  the  anode,  which 
serves  as  the  catalyst.  Therefore,  the  species’  variations  are 
considered  as  mass  fluxes  going  inward  or  outward  from 
the  anode  surface.  Mass  fluxes  of  CO,  C02,  CH4,  H20  on 
anode/fuel  interface  in  a  Av  section  at  position  x  are  thus  in 
the  form  of: 


m 


h2  _ 


=  MH: 


1 


277Tf  Av 


(3-TAx  +  yAx  ^Ax) 


(29) 


mf°  =  M  co 


1 


27ZTf  Av 


C^Ax  jAx) 


(30) 


^°2  =  Mco^ 


jAx 


(31) 


rh^2°  =  M\\2q 


1 


2nrf  Ax 


(^Ax  %Ax  JAx) 


(32) 


m™4  =  Men , - 

*  42jtrfAx 


(  %Ax) 


(33) 


In  a  Av  section  at  position,  x,  the  mass  flux  of  02  on  cath¬ 
ode/air  interface  towards  the  electrolyte  is: 


m 


o2  _ 


=  M() 


1 


27rrair  Av 


(zl2x) 


(34) 


With  the  reacted  mole  numbers  of  CH4  and  CO  determined, 
the  heat  released  from  reforming  and  shift  reaction  is  thus 
in  the  forms  of: 


e^“g  =  A  preforming (35) 

QfZ  =  AHshmyAx  (36) 


2.3.  Flow  and  heat/mass  transfer 


From  the  preceding  discussions,  an  axial- symmetrical 
two-dimensional  (x-r)  computation  domain  is  profiled  as 
shown  in  Fig.  4,  which  includes  two  flow  streams  and  solid 


area  that  locates  the  cell  tube  and  air-inducing  tube.  The 
following  governing  equations,  Eqs.  (37)-(40),  for  momen¬ 
tum,  energy  and  species  transport  are  applicable  universally 
to  the  entire  computation  domain,  however,  zero  velocities 
need  to  be  assigned  to  solid  area  in  the  numerical  treatment. 
As  the  species  fractions  vary  in  the  flow  field,  all  thermal 
and  transport  properties  in  these  equations  are  localized  by 
taking  into  consideration  of  the  local  species’  concentration, 
temperature  and  pressure  [16,17].  This  is  a  specific  point 
for  the  flow,  heat  and  mass  transfer  governing  equations  in 
fuel  cells. 


d(pu)  1  d(rpv) 

dx  r  dr 


(37) 


d(puu)  1  3  (rpvu) 
dx  r  dr 
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d(puv)  1  3  (rpvv) 
dx  r  dr 
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dx  r  dr 


3 
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d(puYf)  1  d{rpvY\) 

dx  r  dr 
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dx 


pDi, 


dYi 

dx 
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1  3 

r  dr 


rpDi, 


dYi 

dr 


m 


(41) 


A  mesh  system  of  at  least  66  in  the  r-direction  and  602  in  the 
v-direction  is  deployed  for  the  computation  domain.  Then, 
all  the  above  equations  are  discretized  by  using  the  finite 
volume  approach  and  the  SIMPFE  algorithm  was  adopted  to 
treat  the  coupling  of  the  velocity  and  pressure  fields  [18,19]. 

The  energy  equation  applied  to  the  solid  components  of 
the  cell  reduces  to  a  heat  conduction  equation  since  zero  ve¬ 
locity  is  assigned  there.  Heat  generation  rates  are  introduced 
into  the  source  term  of  Eq.  (40).  The  heat  conductivity  for 
cell  materials  is  given  in  Table  1 .  The  temperature  between 
the  cell  tube  and  the  air-inducing  tube  might  be  high  enough 
and,  radiation  heat  exchange  might  exist.  Numerical  treat¬ 
ment  for  this  radiation  heat  exchange  is  based  on  the  method 
introduced  in  the  literature  [20]. 

The  boundary  conditions  for  the  momentum,  heat  and 
mass  conservation  equations  are  as  follows.  On  the  symmet¬ 
rical  axis,  or  at  r  =  0:  v  =  0  and  30/3 r  =  0,  where  <fi  repre¬ 
sents  general  variables  except  v.  At  the  outmost  boundary  of 
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Fig.  4.  Computational  domain  pertain  to  one  SOFC  outlined  from  a  cell  stack. 


r  =  /%:  thermally  adiabatic  conditions,  impermeability  for 
species  and  non-chemical  reaction  were  assumed,  that  gives, 
v  =  0  and  30/ dr  =  0,  where  0  represents  general  variables 
except  v.  At  x  =  0:  the  fuel  inlet  has  prescribed  uniform  ve¬ 
locity,  temperature  and  species  mass  fraction,  and  the  solid 
part  has  u  =  0,  v  =  0,  dT/dx  =  0  and  dY\/dx  =  0.  At 
x  =  L  \  the  air  inlet  has  prescribed  uniform  velocity,  temper¬ 
ature  and  species  mass  fraction;  the  gas  exit  part  has  v  =  0, 
du/dx  =  0,  dT/dx  =  0  and  dYi/dx  =  0,  and  all  the  tube  end 
solid  part  has  u  =  0,  v  =  0,  dT/dx  =  0  and  dY\/dx  =  0. 

At  the  interfaces  of  air/support  tube,  r  =  ra ir,  and 
fuel/anode,  r  =  rf,  u  =  0  is  assumed.  However,  in  the 
fuel  flow  passage,  the  mass  flow  rate  increases  along  the 
x-direction  due  to  the  transferring-in  of  oxide  ions.  Simi¬ 
larly,  reduction  of  air  flow  rate  occurs  in  the  air  flow  pas¬ 
sage,  due  to  the  ionizing  of  oxygen  and  ion-transference  to 
the  fuel  side.  Therefore,  a  radial  velocity  exists  at  r  =  rair 
and  r  =  rf 


Vf  = 

Fair 


fuel,  species 

’A 


nfuel 

Px 


air,  species 


P 


air 


r=r  f  ? 


r=ra  ir 


(42) 


The  mass  fraction  at  the  boundaries  of  r  =  rair  and  r  = 
rf  of  all  participating  chemical  components  are  calculated 
with  consideration  of  both  diffusion  and  convection  effects 
[21]: 


m'f  = 


DfairP 


(43) 


«4fUel  =  Hi,  fuel  P 


l  ucl  A 
*  dr 


+  pflYm 


(44) 


It  is  worth  mentioning  that  the  mass/mole  fraction  obtained 
from  this  calculation  will  be  used  for  the  determination  of 
the  partial  pressures  and  thus  the  local  electromotive  forces 
with  Eq.  (4). 


2.4.  Profile  of  the  overall  simulation  scheme  for  a  tubular 
SOFC 


The  species’  mole  compositions,  temperature  and  flow 
rates  of  fuel  and  air  are  the  basic  prescribed  conditions  for 
a  fuel  cell  operation.  In  addition,  the  cell  terminal  voltage 
is  another  decisive  condition  to  designate  a  fuel  cell  operat¬ 
ing  condition.  Other  resultant  parameters  for  a  fuel  cell  will 
follow  these  conditions.  In  the  same  style,  a  numerical  sim¬ 
ulation  can  also  follow  these  prescribed  conditions  to  pre¬ 
dict  the  output  current  and  all  the  other  internal  details,  for 
example,  the  local  temperature,  mass  fractions,  percentage 
of  utilized  fuel,  oxygen,  etc.  However,  instead  of  the  cell 
terminal  voltage,  one  can  also  prescribe  the  output  current, 
together  with  the  basic  prescribed  conditions,  to  designate  a 
fuel  cell  operating  condition.  In  such  a  case,  following  the 
prescribed  conditions,  the  cell  terminal  voltage  and  other  de¬ 
tails  will  be  predicted.  This  is  the  scheme  used  in  this  study. 

It  is  quite  common  in  practice  that  the  output  current  is 
prescribed  in  terms  of  the  average  current  density  of  the  fuel 
cell.  Also,  instead  of  the  flow  rates  of  fuel  and  air,  the  stoi¬ 
chiometric  data  are  prescribed  in  terms  of  the  utilization  per¬ 
centage  of  hydrogen  and  oxygen.  This  kind  of  designation 
of  the  operating  conditions  gives  the  convenience  for  com¬ 
paring  the  fuel  cell  performance  based  on  the  same  level  of 
average  current  density  and  the  hydrogen  and  oxygen  uti¬ 
lization  percentage.  The  inlet  velocities  of  fuel  and  air  are, 
then,  obtainable  in  the  forms  of: 

^cell-^cell  \  RTf 
2F£/H2AH2Afuel/  Pf 


Wfuel  — 


^air  — 


^cell^cell 


RT 


air 


4FUq?  Xq?  Aa  ir 


air 


(46) 


where  Afuei  and  Aair  are  the  cross-sectional  inlet  areas  of 
the  fuel  and  air;  Pf ,  Pa\r  and  7f,  Ta\r  are  the  inlet  pressures 
and  temperature  of  the  fuel  and  air  flows,  respectively;  Xh2 
and  Xo2are  the  mole  fractions  of  hydrogen  in  the  fuel  and 
oxygen  in  the  air,  respectively. 
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The  computation  process  is  highly  iterative  and  interdigi- 
tated  in  nature.  As  the  first  step,  the  local  temperature,  pres¬ 
sure,  species’  mass  fraction  are  used  in  the  network  circuit 
analysis  to  obtain  the  cell  terminal  voltage,  local  current 
across  the  electrolyte  and  thus  the  local  species’  transfer 
fluxes  and  heating  rates.  In  the  second  step,  the  local  tem¬ 
perature,  pressure  and  species’  mass  fractions  are,  in  turn, 
obtained  through  solution  of  the  governing  equations  un¬ 
der  new  internal  and  boundary  conditions  determined  by  the 
latest- available  species’  flux  and  heating  rate.  The  two  steps 
iterate  interactively  until  convergence  is  obtained. 

3.  Results  and  discussion  for  validation  computation 

Experimental  data  and  operating  conditions  for  tubular 
SOFCs  from  different  papers,  respectively,  by  Hagiwara 
et  al.  [22],  Hirano  et  al.  [7],  Singhal  [23,24]  and  Tomlins 
and  Jaszar  [25]  have  been  referred  for  the  validation  of  the 
present  numerical  model.  The  fuel  used  by  Hirano  et  al.  [7] 
had  components  of  H2,  H2O,  CO  and  CO2.  Thus,  there  is 
a  water-shift  of  the  carbon  monoxide  in  the  fuel  cell  to  be 
considered.  The  fuel  used  by  other  researchers  [22-25]  had 
components  of  H2  and  H2O  and  there  is  only  the  electro¬ 
chemical  reaction  to  be  considered.  The  dimensions  of  the 
three  different  solid  oxide  fuel  cells  tested  in  their  studies 
are  listed  in  Table  2,  in  which  the  mesh  size  adopted  in 
our  numerical  computation  is  also  given.  Correspondingly, 
Table  3  lists  the  species  mole  fraction  and  temperature  of 
fuel  and  air  in  their  test,  which  are  the  prescribed  conditions 
in  this  computation.  The  experimental  data  for  these  SOFCs 
in  different  structural  dimensions  and  operating  conditions 


0  100  200  300  400  500  600 

Current  Density  (  mA/cm2) 

Fig.  5.  Prediction  and  testing  results  of  cell  voltage  vs.  current  density. 
(Operating  pressure  of  the  cell  tested  by  Hagiwara  et  al.  [22]  and  Hirano 
et  al.  [7]  are  1.013  x  105  Pa,  and  that  by  Tomlins  and  Jaszar  [25]  is 
5.065  x  105  Pa). 

is  expected  to  facilitate  a  wide  benchmark  for  validation  of 
the  numerical  model  presented  in  this  study. 

3.1.  The  SOFC  terminal  voltage  versus  current  density  for 
different  cells 

Fig.  5  shows  the  numerically  predicted  cell  terminal  volt¬ 
age  and  the  experimental  data  in  the  literature  at  different 
cell  current  densities.  The  relative  deviation  between  the 


Table  2 

Dimensions  of  the  cited  SOFCs  and  mesh  size  used  in  this  computation 


Data  sequence:  outer  diameter  (mm)/thickness  (mm)/length  (mm) 


[22] 

[7] 

[23-25] 

Air-inducing  tube 

7.00/1.00/485 

6.00/1.00/290 

12.00/1.00/1450 

Support  tube 

— 

13.00/1.50/300 

— 

Cathode 

15.72/2.20/500 

14.40/0.70/300 

21.72/2.20/1500 

Electrolyte 

15.80/0.04/500 

14.48/0.04/300 

21.80/0.04/1500 

Anode 

16.00/0.10/500 

14.68/0.10/300 

22.00/0.10/1500 

Fuel  boundary 

18.10/-/500 

16.61/-/300 

24.87/-/1500 

Grid  number  (r  x  x) 

66  x  602 

66  x  602 

66  x  1602 

Table  3 


Species’  mole  fraction,  utilization  percentage  and  temperature 


Air  [02  (%)-£/o2/N2  (%)/T  (°C)] 

Fuel  [H2  (%)-C/h2/H20  (%)/CH4  (%)/CO  (%)/C02  (%)/T  (°C)] 

[22] 

21.00-17.00/79.00/600.0 

98.64-85.00/1.36/0/0/0/900.0 

[7] 

2 1 .00-25.00/79.00/600.0“ 
21.00-25.00/79.00/400.0b 

55.70-80.00/27.70/0/10.80/5.80/800.0 

55.70-80.00/27.70/0/10.80/5.80/800.0 

[24,25] 

21.00-17.00/79.00/600.0 

98.64-85.00/1.36/0/0/0/800.0 

a  Case  of  current  density  =185  mA/cm2 . 
b  Case  of  current  density  =  370  mA/cm2. 
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model-predicted  cell  voltages  and  the  experimental  data  is  no 
larger  than  1.0%  for  the  cell  tested  by  Hirano  et  al.  [7],  5.6% 
for  that  by  Hagiwara  et  al.  [22],  and  6.0%  for  that  by  Tomlins 
and  Jaszar  [25].  Such  a  good  agreement  for  the  terminal 
voltages  between  the  model-prediction  and  the  experiment 
from  different  researchers  is  a  strong  demonstration  that  the 
present  numerical  model  is  quite  reliable. 

A  comparison  of  the  fuel  cell  performance  in  terms  of  the 
cell  voltage  versus  the  cell  current  density  for  the  simulated 
results  of  the  three  SOFCs  in  different  diameter  and  length 
may  be  necessary.  It  is  found  from  Fig.  5  that,  under  the  same 
cell  current  density,  the  cell  voltage  of  the  SOFC  tested  by 
Hagiwara  et  al.  [22]  is  the  highest  and  that  by  Hirano  et  al. 
[7]  is  the  lowest.  The  mole  fraction  of  hydrogen  is  low  in  the 
fuel  for  the  SOFC  tested  by  Hirano  et  al.  [7],  which  might  be 
the  major  reason  that  it  has  the  lowest  cell  voltage.  Because 
the  current  must  be  collected  circumferentially  in  a  tubular 
type  fuel  cell,  the  large  diameter  of  cell  tube  by  Singhal  [24] 
and  Tomlins  and  Jaszar  [25]  will  lead  to  a  longer  current 
pathway  and  thus  their  cell  voltage  is  lower  than  that  by 
Hagiwara  et  al.  [22]  even  though  the  data  shown  in  Fig.  5 
was  at  a  pressurized  operation  of  5.065  x  105  Pa,  which,  in 
fact,  helps  to  improve  the  cell  voltage. 

Under  a  constant  current  density  of  300mA/cm2,  the  cell 
voltage  and  power  increase  with  the  increasing  operating 
pressure,  as  seen  in  Fig.  6.  The  agreement  between  our  simu¬ 
lation  results  and  the  experimental  data  by  Singhal  [24]  is 


(b)  Operating  Pressure  (atm) 


Fig.  6.  Prediction  and  testing  results  of  the  cell  voltage  and  power  vs. 
operating  pressure.  (Experimental  data  was  from  Singhal  [24]  under  cell 
current  density  of  300m A/cm2.) 


Fig.  7.  Cell  temperature  distribution  from  prediction  and  experimental 
testing.  (Experimental  data  are  from  Hirano  et  al.  [7]  at  cell  current 
density  of  185mA/cm2.) 


quite  good,  showing  a  maximum  deviation  of  7.4%  at  low 
operating  pressures.  When  the  operating  pressure  increases 
from  1.013  x  105  to  5.065  x  105  Pa,  the  cell  output  power 
shows  a  significant  improvement  of  9%.  However,  raising 
the  operating  pressure  becomes  less  effective  for  improv¬ 
ing  the  output  power  when  the  operating  pressure  is  high, 
for  example,  the  cell  output  power  has  only  a  6%  increase 
when  the  operating  pressure  increases  from  5.065  x  105  to 
1.520  x  106  Pa.  This  is  because  that  the  operating  pressure 
contributes  to  the  cell  voltage  in  a  logarithmic  manner.  Nev¬ 
ertheless,  the  cell  output  power  can  be  improved  significantly 
by  pressurized  operation.  Increasing  the  operating  pressure 
from  1.013  x  105  to  1.520  x  106Pa,  the  cell  output  power 
can  have  an  increment  of  15.8%. 

3.2.  Results  of  the  cell  temperature  distribution 

Because  the  measurement  of  temperature  in  a  SOFC  is 
very  difficult,  only  three  experimental  data,  the  temperature 
at  two  ends  and  in  the  middle  of  the  cell  tube,  was  available 
from  the  work  by  Hirano  at  al.  [7].  Shown  in  Fig.  7  is  the 
simulated  cell  temperature  distribution  for  the  SOFC  tested 
by  Hirano  et  al.  [7].  The  agreement  of  the  simulated  data 
and  the  experimental  results  is  rather  good  in  the  middle, 
where  the  hotspot  locates,  and  relatively  larger  deviation 
appears  at  the  two  ends  of  the  cell.  For  further  compari¬ 
son,  more  experimental  data  tested  for  different  SOFCs  is 
expected  in  the  future. 

The  predicted  temperature  distribution  for  the  fuel  cells 
tested  by  Hagiwara  et  al.  [22],  and  Tomlins  and  Jaszar  [25]  is 
given  in  Fig.  8.  Generally,  the  two  ends  of  the  cell  tube  have 
a  lower  temperature  than  the  middle  of  the  cell  tube.  How¬ 
ever,  at  low  current  densities,  the  hotspot  is  located  closer  to 
the  closed  end  of  the  cell.  With  increasing  current  density, 
the  hotspot  shifts  to  the  open-end  side  and  its  temperature 
also  decreases,  which  makes  the  uniformity  of  temperature 
distribution  along  the  fuel  cell  better.  The  heat  transfer  be¬ 
tween  the  cooling  air  and  the  cell  tube  at  the  cell  closed-end 
region  is  like  a  jet  impingement  heat  transfer,  although  the 
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Fig.  8.  Predicted  cell  tube  temperature  distribution  under  different  current 
density:  (a)  the  cell  tested  by  Hagiwara  et  al.  [22];  (b)  the  cell  tested  by 
Singhal  [24]  and  Tomlins  and  Jaszar  [25]. 


exit  velocity  from  the  air-inducing  tube  is  quite  low.  There¬ 
fore,  the  velocity  of  the  exit  air  from  the  air-inducing  tube 
affects  the  heat  transfer  coefficient  significantly.  For  the  high 
current  density  case,  the  air  flow  rate  becomes  larger  ac¬ 
cordingly,  thus  the  heat  transfer  coefficient  between  the  air 
and  the  fuel  cell  closed-end  region  is  increased.  This  can 
prohibit  the  temperature  level  of  the  closed-end  region  of 
the  fuel  cell  significantly.  Since  the  air  receives  much  heat 
at  the  closed-end  region,  its  cooling  for  the  fuel  cell  in  the 
downstream  region  becomes  weak  and  thus,  the  uniformity 
of  cell  temperature  distribution  becomes  much  better  when 
at  high  current  densities. 

3.3.  Details  of  flow,  temperature  and  concentration  fields 

Fig.  9  shows  the  flow  and  temperature  fields  for  the 
SOFC  tested  by  Hirano  et  al.  [7]  at  a  current  density  of 


185mA/cm2.  The  air  speed  in  the  air-inducing  tube  has  a 
slight  acceleration  because  the  air  absorbs  heat  and  expands 
in  this  flow  passage.  After  leaving  the  air-inducing  tube,  the 
air  gives  an  impingement  to  the  closed  end  of  the  fuel  cell 
and  then  flows  backward  to  the  outside.  In  this  pathway,  the 
air  obtains  heat  from  the  heat-generating  fuel  cell  tube  and 
transfers  this  heat  to  the  cold  air  in  the  air-feeding  tube.  It 
is  easy  to  understand  that  the  electrochemical  reaction  at 
the  closed  end  of  the  fuel  cell  is  strong  because  the  concen¬ 
trations  of  fuel  and  air  are  both  high  there.  Therefore,  the 
heat  generation  due  to  Joule  heating  and  entropy  change  of 
the  electrochemical  reaction  is  higher  at  the  upstream  of  the 
fuel  path.  However,  the  fuel  cell  closed  end  does  not  have 
the  highest  temperature,  therefore,  it  is  believed  that  the 
cooling  of  air  to  the  closed  end  of  the  fuel  cell  is  responsi¬ 
ble  for  this.  After  being  heated  at  the  closed-end  region,  air 
temperature  is  high  and  its  cooling  ability  for  the  cell  tube 
is  low  when  it  is  in  the  annulus  between  the  air-inducing 
tube  and  the  cell  tube.  At  the  cell  open-end  region,  the  air  in 
the  annulus  can  transfer  heat  to  the  incoming  cold  fresh  air 
in  the  air-inducing  tube  and  this  will  help  it  to  cool  the  fuel 
cell  tube.  From  this  air  flow  arrangement,  the  hotspot  tem¬ 
perature  of  the  cell  tube  might  mostly  occur  in  the  center  of 
the  longitudinal  direction  of  the  cell  tube.  The  air  flow  has 
two  passes,  incoming  in  the  air-inducing  tube  and  outgoing 
in  the  annulus  between  the  air-inducing  tube  and  cell  tube, 
therefore,  the  heat  exchange  in  between  the  two  passes  al¬ 
lows  the  air  to  mitigate  its  temperature  fluctuation  along  the 
whole  air  path  and  thus  the  temperature  field  in  the  fuel  cell 
might  be  maintained  relatively  uniform.  Nevertheless,  the 
heat  generation,  air  and  fuel  temperature  and  air-cooling  to 
the  fuel  cell  will  collectively  affect  the  temperature  field  in 
the  fuel  cell.  Thus,  the  hotspot  position  in  a  cell  tube  might 
shift  more  or  less  away  from  the  center  region  depending 
on  the  operating  condition  of  the  fuel  cell. 

Fig.  10  shows  the  gas  species’  molar  fraction  contours  in 
even-difference  for  the  same  SOFC  under  the  same  operat¬ 
ing  condition  as  discussed  in  Fig.  9.  In  the  air  path,  oxygen 
consumption  at  the  cell  closed-end  region  is  relatively  large, 
which  leads  to  more  dense  contour  lines.  The  contour  shape 
of  oxygen  also  indicates  a  relatively  large  difference  of  the 
molar  fraction  between  the  bulk  flow  and  the  wall  of  cath¬ 
ode/air  interface.  This  also  implies  that  the  mass  transport 
resistance  on  airside  might  be  dominant  in  lowering  the  cell 
performance.  Feeding  more  air  than  needed  is  the  way  al¬ 
ready  well  accepted  to  in  fuel  cell  technology. 

The  consumption  by  electrochemical  reaction  and  gener¬ 
ation  from  water- shift  of  CO  collectively  determines  the  hy¬ 
drogen  budget.  Since  the  consumption  takes  the  upper  hand, 
the  hydrogen  molar  fraction  decreases  along  the  fuel  stream. 
Corresponding  to  this  situation  for  hydrogen,  the  consump¬ 
tion  due  to  water-shift  and  production  due  to  electrochemical 
reaction  makes  the  water  vapor  increase  gradually  along  the 
fuel  stream.  The  water-shift  of  CO  proceeds  gradually  along 
the  fuel  path,  and  thus  the  mole  fraction  of  CO  decreases 
and  the  CO2  increases.  The  shape  of  the  contour  lines  of  the 
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Fig.  9.  Simulation  results  of  the  flow  and  temperature  fields  for  the  cell  tested  by  Hirano  et  al.  [7].  (Cell  current  density  is  185mA/cm2.) 
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Fig.  10.  Simulation  results  of  the  species’  molar  fraction  contours  for  the  cell  tested  by  Hirano  et  al.  [7].  (Cell  current  density  is  185mA/cm2.) 
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Fig.  11.  Simulation  results  of  species’  molar  flow  rate  varying  in  the  fuel 
channel  for  the  cell  tested  by  Hirano  et  al.  [7].  (Cell  current  density  is 
185  mA/cm2.) 

species  in  the  fuel  path  is  relatively  flat  from  the  cell  wall 
to  the  bulk  flow.  This  is  an  indication  that  mass  diffusion  in 
the  fuel  channel  is  relatively  strong.  For  further  understand¬ 
ing  of  the  variation  of  the  gas  species,  Fig.  1 1  shows  their 
molar  flow  rate  variation  along  the  fuel  path.  In  one  third 
of  the  length  from  fuel  inlet,  the  hydrogen  flow  rate  shows 
a  faster  decrease  and  the  water  flow  rate  shows  a  faster  in¬ 
crease,  indicating  a  stronger  reaction  in  the  upstream.  The 
flow  rate  of  CO  and  CO2  vary  roughly  in  a  linear  style  and 
a  small  amount  of  CO  still  exists  in  the  waste  gas. 

4.  Conclusion 

The  present  model  assumed  that  the  heat  and  mass  species 
do  not  exchange  in  between  one  cell  and  its  neighboring 
cells  when  the  tubular  SOFCs  is  working  in  a  cell  stack. 
A  surrounding  fuel  flow  space  controllable  by  the  inter¬ 
ested  cell  was  then  outlined,  for  which  zero  heat  and  mass 
fluxes  were  assumed  at  the  boundary  in  the  neighborhood 
with  other  cells.  A  two-dimensional  numerical  model  for 
the  flow,  heat/mass  transfer  and  chemical/electrochemical 
performance  in  such  a  cell  operating  in  a  stack  was  devel¬ 
oped.  Numerical  simulation  was  conducted  for  three  dif¬ 
ferent  tubular  SOFCs  having  practical  operational  results 
published  by  different  researchers.  The  model-predicted  cell 
voltage  agrees  with  the  experimental  data  with  a  deviation 
of  less  than  7.4%  for  all  the  experimental  data.  The  fuel  cell 
temperature  from  the  prediction  also  agrees  with  the  test 
data  satisfactorily. 

It  was  found  that  the  two  ends  of  fuel  cell  have  a  rel¬ 
atively  low  temperature.  The  hotspot  locates  closer  to  the 
closed  end  of  the  cell  tube  when  the  operating  current  den¬ 
sity  is  low.  With  increasing  current  density,  the  hotspot  shifts 
to  the  open-end  side  of  the  cell  and  also  the  hotspot  tem¬ 


perature  decreases  which  makes  the  uniformity  of  temper¬ 
ature  distribution  better  along  the  fuel  cell  tube.  The  molar 
fraction  contour  of  oxygen  in  the  air  channel  shows  a  large 
difference  in  the  oxygen  concentration  between  the  bulk  air 
flow  and  the  cathode/air  interface.  In  the  fuel  channel,  the 
species  molar  fraction  contour  lines  are  rather  flat,  showing 
a  relatively  smaller  difference  between  bulk  flow  and  the 
anode/fuel  interface. 

References 

[1]  S.  Srinivasan,  R.  Mosdale,  R  Stevens,  C.  Yang,  Annu.  Rev.  Energy 
Environ.  24  (1999)  281-328. 

[2]  F.J.  Gardner,  Proc.  Inst.  Mech.  Eng.  A  211  (1997)  367-379. 

[3]  S.P.  Harvey,  H.J.  Richter,  ASME  J.  Energy  Resour.  Technol.  116 
(1994)  312-318. 

[4]  K.  Suzuki,  H.  Iwai,  J.H.  Kim,  P.-W.  Li,  K.  Teshima,  Proceedings  of 
the  12th  International  Heat  Transfer  Conference,  Grenoble,  18-23 
August  2002,  vol.  1,  2002,  pp.  403-414. 

[5]  E.F.  Sverdrup,  CJ.  Warde,  R.L.  Eback,  Energy  Convers.  13  (1973) 
129-136. 

[6]  N.F.  Bessette,  W.J.  Wepfer,  Trans.  ASME  J.  Energy  Resour.  Technol. 
117  (1995)  43-49. 

[7]  A.  Hirano,  M.  Suzuki,  M.  Ipponmatsu,  J.  Electrochem.  Soc.  139 
(1992)  2744-2751. 

[8]  P.  Aguiar,  D.  Chadwick,  L.  Kershenbaum,  Chem.  Eng.  Sci.  57  (2002) 
1665-1677. 

[9]  J.R.  Ferguson,  J.M.  Fiard,  R.  Herbin,  J.  Power  Sources  58  (1996) 
109-122. 

[10]  W.M.  Kays,  M.E.  Crawford,  Convective  Heat  and  Mass  Transfer, 
third  ed.,  McGraw-Hill,  New  York,  1993. 

[11]  P.-W.  Li,  L.  Schaefer,  Q.M.  Wang,  M.K.  Chyu,  in:  Paper  Presented 
at  the  ASME  International  Mechanical  Engineering  Congress  and 
Exposition,  Paper  No.  IMECE2002-32564,  17-22  November  2002, 
New  Orleans,  LA,  USA. 

[12]  E.A.  Liese,  R.S.  Gemmen,  F.  Jabbari,  J.  Brouwer,  in:  Proceedings 
of  the  International  Gas  Turbine  and  Aero-engine  Congress  and 
Exhibition,  Indianapolis,  IN,  USA,  ASME  Paper  99-GT-360,  1999. 

[13]  S.H.  Chan,  K.A.  Khor,  Z.T.  Xia,  J.  Power  Sources  93  (2001)  1 30— 
140. 

[14]  S.  Ahmed,  C.  McPheeters,  R.  Kumar,  J.  Electrochem.  Soc.  138  (9) 
(1991)  2712-2718. 

[15]  A.F.  Massardo,  F.  Lubelli,  Paper  Presented  at  the  International 
Gas  Turbine  and  Aero-engine  Congress  and  Exhibition,  Paper  No. 
98-GT-577,  1998. 

[16]  R.H.  Perry,  D.W.  Green,  Perry’s  Chemical  Engineer’s  Handbook, 
seventh  ed.,  McGraw-Hill,  New  York,  1986. 

[17]  P.F.  Incropera,  D.P  DeWitt,  Introduction  to  Heat  Transfer,  third  ed., 
Wiley,  New  York,  1996. 

[18]  S.V.  Patankar,  Numerical  Heat  Transfer  and  Fluid  Flow, 
McGraw-Hill,  New  York,  1980. 

[19]  K.C.  Karki,  S.V.  Patanka,  AIAA  J.  27  (9)  (1989)  1167-1174. 

[20]  C.  Beckermann,  T.F.  Smith,  J.  Numer.  Heat  Transfer,  Part  B  23 
(1993)  127-133. 

[21]  E.R.G.  Eckert,  R.M.  Drake,  Heat  and  Mass  Transfer,  McGraw-Hill, 
New  York,  1959. 

[22]  A.  Hagiwara,  H.  Michibata,  A.  Kimura,  M.P.  Jaszcar,  G.W.  Tomlins, 
S.E.  Veyo,  in:  Proceedings  of  the  Third  International  Fuel  Cell 
Conference,  D2-4,  1999,  pp.  365-368. 

[23]  S.C.  Singhal,  Electrochem.  Soc.  Proc.  19-99  (1999)  39-51. 

[24]  S.C.  Singhal,  Solid  State  Ionics  135  (2000)  305-313. 

[25]  G.W.  Tomlins,  M.P  Jaszar,  in:  Proceedings  of  the  Third  International 
Fuel  Cell  Conference,  D2-4,  1999,  pp.  369-372. 


